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We analyze the limits on resonant particle production dur- 
ing inflation based upon the power spectrum of fluctuations in 
matter and the cosmic microwave background. We show that 
such a model is consistent with features observed in the mat- 
ter power spectrum deduced from galaxy surveys and damped 
Lyman-Qf systems at high redshift. It also provides an alter- 
native explanation for the excess power observed in the power 
spectrum of the cosmic microwave background fluctuations in 
the range of 1000 < I < 3500. For our best-fit models, epochs 
of resonant particle creation reenter the horizon at wave num- 
bers of kt. ~ 0.4 and/or 0.2 (ft Mpc~^). The amplitude and 
location of these features correspond to the creation of fermion 
species of mass '-^ 1 — 2 Mp; during inflation with a coupling 
constant between the inflaton field and the created fermion 
species of near unity. Although the evidence is marginal, if 
this interpretation is correct, this could be one of the first 
observational hints of new physics at the Planck scale. 



I. INTRODUCTION 

Analysis of the power spectrum of fluctuations in the 
the large-scale distribution of matter (cf. [2-4]), together 
with fluctuations in the cosmic microwave background 
(CMB) (cf. [1]) provides powerful constraints on the 
physics of the very early universe. The most popular 
account for the origin of both power spectra is based 
upon quantum fluctuations generated during the infla- 
tionary epoch [5]. Subsequently, acoustic oscillations of 
the photon-baryon fluid distort this to produce the ob- 
served features in the angular power spectrum of tem- 
perature fluctuations in the CMB and the spatial power 
spectrum of matter density fluctuations. The two power 
spectra are different in that the cosmic microwave back- 
ground is sensitive to the baryonic material, while the 
matter power spectrum probes the dominating dark mat- 
ter. 

In this context, there now exist determinations of the 
matter power spectrum on small angular scales due to re- 
cent SDSS [3,17], and 2dF [18] galaxy surveys as well as 
analysis of the Lyman-a forest [19,20] at higher redshift 
and even smaller scales. This latter determination is par- 
ticularly facilitated by the fact that at high redshift, the 
Lyman-o; absorption systems are still within the quasi- 
linear regime so that an inference of the primordial power 
spectrum is relatively straightforward. The overall mat- 



ter power spectrum deduced in this way is consistent with 
a standard ACDM cosmology, but as shown in Figure 1 
there is marginal evidence of a peculiar feature beginning 
near k ^ 0.6 h Mpc^^ which is not easily explained away 
by systematic errors [20]. There is also at least a pos- 
sibility for structure to exist in the region near k ~ 0.3 
h Mpc~^. In this paper we consider that such features 
may be a part of the primordial spectrum generated by 
new physics near the end of the inflation epoch. 
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FIG. 1. Comparison of the observed galaxy cluster func- 
tion from SDSS [17], 2dF [18], and Lyman-a [19] with the 
spectrum implied from the fits to the matter power spectrum 
with (solid line) and without (dashed line) resonant particle 
creation during inflation as described in the text. 

Regarding the CMB, data [6] from WMAPhave placed 
stringent constraints on the cosmological parameters rel- 
evant to the observed power spectrum in the range of 
multipoles up to Z < 1000. The simplest prediction of 
the standard ACDM cosmology is that the primordial 
power spectrum should be strongly damped at higher 
multipoles due to photon diffusion (Silk damping) and 
because more than one perturbation can fit within the 
depth of the surface of last scattering, and hence, be 
washed out. 

Of interest for the present work, however, is the recent 
accumulation of observations of the CMB power spec- 
trum in the range 1000 < I < 7000 by various groups 
{CBI [7-9], ACBAR [10], BIMA [11], and VSA [12]). 
These observations on the smallest angular scales are of 



particular interest as a test of this basic prediction of 
the standard inflation/photon-decouphng paradigm. The 
current data are summarized in Figure 2 from which it is 
clear that the deduced power spectrum increases rather 
than decreases for large multipoles, particularly in the 
range 2000 < I < 3500. 
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FIG. 2. Upper figure shows the CMB WMAP, ACBAR, 
and CBI data in the range / = 10 — 5000. The dashed line 
is the CMB power spectrum computed using the standard 
WMAP cosmological parameters. The thick solid line is for 
a best-fit to these data for a model with resonant particle 
production included. The lower figure shows an expanded 
view of the ACBAR, CBI, VSA, and BIMA data in the range 
of I = 1000 - 7000. The dashed line shows the standard 
WMAP result without the SZ effect included. The dot-dashed 
line shows our fit SZ contribution (from the analytic halo 
model of Ref. [14]). The solid line shows the best fit to these 
data sets with resonant particle production and the SZ effect 
included. The dot-dot-dashed line shows the resonant particle 
creation component without the SZ contribution. 

The most likely interpretation (cf. [13,14]) of an ex- 
cess power at high multipoles is a manifestation of the 
scattering of CMB photons by hot electrons in clusters 
known as the thermal Sunyaev-Zeldovich (SZ) effect [15]. 
Although some contribution from the SZ effect undoubt- 



edly occurs, it is not yet conclusively established that 
this is the only possible interpretation (cf. [16]) of ex- 
cess power on small angular scales. Indeed, it has been 
deduced [13,14] that explaining the excess power in the 
observed CMB spectrum requires that the mass density 
fluctuation amplitude parameter, a^, be near the upper 
end of the range of the values deduced via other indepen- 
dent means. Another significant effect is the large cosmic 
variance. The uncertainties on the SZ contribution are 
highly non-Gaussian. They thus depend strongly upon 
whether or not the observed field happens to contain a 
group or cluster of galaxies. Hence, there is large uncer- 
tainty in any deduced SZ effect. In view of these uncer- 
tainties, an exploration of other possible explanations for 
this power excess seems warranted. 

Indeed, alternative explanations exist for the genera- 
tion of features in both the matter and CMB power spec- 
tra on small angular scales. For example, a flattening of 
the inflation generating effective potential near the end of 
inflation covild produce such distortions on small angular 
scales [21]. 

In this paper, however, we consider another alterna- 
tive originally proposed in [22] whereby such distortions 
might arise from the resonant production of particles dur- 
ing large field inflation. -'^ This interpretation has the in- 
triguing aspect that, if correct, an opportunity emerges 
to use the CMB and matter power spectra as probes of 
the Planck scale {Mpi ~ 10^^ GeV) particle spectrum. 

The prototypical scenario is that there is an inflaton 
(j) which controls the mass of a fermion t/i or a boson x 
through the coupling 



i^in.t — 



-[Mf - Mp,/(^)]# - [Ml - Mlig{^)]x^ 



(1) 



where Mf and Mb are fermion and boson masses, respec- 
tively, which are assumed to be of order Mpi . Since by the 
deflnition of large fleld inflationary scenario 4>/Mpi > 1 
and (j)/Mpi varies at least an order of magnitude, any 
natural order unity functions f{x) and g(x) will lead 
to a cancellation of the masses during the evolution 
of slow roll inflation. Due to this cancellation, when 
\Mf - Mpifl < H and/or \Mb - Mpig\ < H, the effec- 
tive mass of ip and x is varying nonadiabatically (since 
the inflaton field variation time rate is also typically not 
too far from the Hubble expansion rate H).^ During 



^ Large field infiation is a slow roll infiationary scenario in 
which the infiaton scalar field has a classical field value ex- 
ceeding the Planck scale. 

^More precisely, for quadratic monomial inflaton potentials, 
the nonadiabaticity condition ^ > 1 leads approximately to 



this nonadiabatic period, there is efficient particle pro- 
duction which can be thought of as a kind of resonant 
particle production. 

In [22], the effect of the resonant fermionic particle 
production was taken into account neglecting the nona- 
diabatic effects on the modes outside of the horizon. This 
leads to a bump-like structure in the primordial power 
spectrum. However, [23] considered the nonadiabatic ef- 
fects of bosonic particle production on the modes outside 
the horizon. Their claim is that the primordial power 
spectrum is modified by a steplike structure rather than 
a bumplike structure. Given that their analysis was for 
bosonic particle production while we are going to deal 
with fermionic production following [22],'^ and given that 
the subtle effect of [23] (which was only approximately 
computed numerically) requires more detailed attention 
to unambiguously establish the change in the qualitative 
behavior of the primordial spectral feature, we will in this 
paper consider only the semianalytic results of [22] when 
fitting to data. However, because of this omission of the 
nonadiabatic effect, this result should only be considered 
indicative of the type of constraints one can obtain from 
the latest data. A detailed analysis of the nonadiabatic 
effects will be considered in a future publication. 

Planck-scale mass particles generically exist in Planck 
scale compactification schemes of string theory from the 
Kaluza-Klein states, winding modes, and the massive 
(excited) string modes. Hence, the existence of Planck- 
scale mass particles which couple to the inflaton is a 
generic situation. What is perhaps not generic in the 
scenario considered in this paper is that the mass of the 
Planck-scale particle lies in the e-fold range of the infla- 
ton accessible to observation. 

Even with all of these favorable assumptions, we find 
only marginal hints for the existence of such a parti- 
cle from fits to current data. More specifically, we find 
marginal hints for a bump at a scale of fc ~ 0.4 (and 0.2) 
h Mpc~^ in the primordial spectrum due to production of 
fcrmions of mass m w 2Mpi with an 0(1) weak coupling 
to the inflaton. 

The order of presentation will be as follows. In section 
II, we briefly review and clarify the semianalytic results 
of [22] and then set up the parameterization of the pri- 



H > y{Mf 
for fermion mass term where 
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mordial spectrum fitting function. In sections HI and 
IV, we briefly discuss the matter power spectrum param- 
eterization and the CMB power spectrum (including the 
SZ effect which is the more conservative interpretation of 
the excess power at large I values). Section V is devoted 
to the discussion of the fitting procedure and the fit re- 
sults to the matter and the CMB spectrum. Finally, we 
conclude in section VI. 



II. INFLATION RESONANT PARTICLE 
PRODUCTION 

In the basic inflationary picture, a rapid early expan- 
sion of the universe is achieved through the vacuum en- 
ergy from an inflaton field. In the minimal extension from 
the basic picture considered here, the inflaton is postu- 
lated to couple to at least one massive particle whose 
mass is order of the of the inflaton fleld value. This par- 
ticle is then resonantly produced as the field obtains a 
critical value during infiation. If even a small fraction of 
the energy in the inflaton field is extracted in this way, it 
can produce features in the primordial power spectrum. 
In particular, there will be excess power in the spectrum 
at the angular scale corresponding to when the epoch of 
resonant particle creation crossed the Hubble radius. 

In the simplest slow roll approximation [5] for the gen- 
eration of density perturbations during infiation, the am- 
plitude, Snik), of a density fluctuation when it crosses 
the Hubble radius is just. 



Snik) 






(2) 



_ 47v\g-\Mf/M^t)\ fMf 
^ \g'\ \M,i ^ 



where H is the expansion rate, and (j) is the velocity of the 
inflaton field when the comoving wave number k crosses 
the Hubble radius during inflation. If resonant particle 
production drains energy from the inflaton field, then 
the conjugate momentum in the field (j) decreases. This 
causes an increase in Snik) (primordial power spectrum) 
for those wave numbers which exit the horizon during the 
resonant particle production epoch. 

Of course when (j) is changing due to particle produc- 
tion, (f) may not be negligible, resulting in corrections to 
Eq. (2) . In [22] , this correction was considered and found 
to be << 20% ** for the particle production of interest 
for the fits of this paper. 

The infiaton field is then postulated to have a simple 
Yukawa coupling to a fermion field tp of mass m in the 
form, 



Cy = —\4>ilJip 



(3) 



Dealing with fermionic particle production allows us to 
avoid dealing with nonminimal phase transition effects that 
can arise due to the possible nontrivial ground states of the 
bosonic field. 



*This fraction refers to the fraction of the particle production 
effect, not the entire power spectrum amplitude. 



Including this new coupling, the equation of motion for 
the inflaton field becomes 



3H(j) 



dV 



NXixp-iP) = 



(4) 



for N fcrmions of mass m coupled to the inflaton. The 
effective mass of the fermion is M{(1)) ^ m — \(j), which 
vanishes for a critical value of the inflaton field, 0, = 
to/A. Resonant fermion production will then occur in a 
narrow range of inflaton field amplitude around (j) ~ (j),,. 
As in [22] we label the epoch at which particles are 
created by an asterisk. So the cosmic scale factor is la- 
beled a* at the time i» at which resonant particle pro- 
duction occurs. Considering a small interval around this 
epoch, one can treat H = iJ, as approximately constant 
(slow roll inflation). The number density n of parti- 
cles can be taken as zero before t* and afterwards as 
n — n^[ai, / a{t)Y . The fermion vacuum expectation value 
can thus be written. 



(V^V) 



« n* 6l(t-t*)exp[-3iJ4f 



t*)] 



(5) 



Now inserting this relation into the equation of motion 
for the inflaton field (Eq. (4)), one can obtain the change 
in the inflaton field evolution (p due to particle creation. 



(t>{t > U) ^ (i>{t > U)x^o 

+ iVAn*6l(t-t*)exp[-3iJ*(f-i*)] 



(6) 



Inserting this into Eq. (2), a very good analytic ap- 
proximation to the effect of the particle creation on the 
perturbation spectrum can be obtained [22], 
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(7) 



The scale factor a relates [24] to the physical wave num- 
ber fc by, 
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where ao-ffo ~ {h/3000) Mpc~i denotes the present co- 
moving Hubble scale. The subscript 'fc' indicates the in- 
flaton effective potential value when a particular wave 
number k crosses the Hubble radius during inflation 
(k = aH). The quantities aend and Vend are the scale 
factor and effective inflaton potential at the end of infla- 
tion, and proh is the matter energy density after reheating 
to the standard hot big bang Friedmann cosmology. This 
expression assumes that instantaneous transitions occur 



between the various regimes, and that the universe be- 
haves as if matter-dominated during reheating. 

Using this relation between scale factor and k, the per- 
turbation spectrum (Eq. (7)) can be reduced [22] to a 
simple two-parameter function. 



SH{k) = 



[SHJk)]^^, 



l-9{k- k^)A{kJkf ln(fc/fc,) 



(9) 



where the coefficient A and characteristic wave number 
A:* (k/kf: > 1) can be fit to the observed power spectra. 
(Note that this A is different from the A coefficient of 
[22].) 

The values of A and fc* determined from observation di- 
rectly relate to the inflaton coupling A and fermion mass 
TO,, for a given inflation model. When the back reaction 
is not important, we can write 



A 



'^NXn^H-^ 



(10) 



which uses the approximation [22,25-27] for the particle 
production Bogoliubov coefficient to be 



\Pk\'^ = exp 



-TTk^ 

a*A|(/), 



(11) 



Note that this approximation does not depend on the 
particular form of the inflaton potential. Its main as- 
sumptions are only that the particle mass being pro- 
duced is nearly negligible during the resonant produc- 
tion (which is an excellent assumption for our case) and 
that the effective-mass time variation due to the infla- 
ton time variation dominates over the contribution due 
to the FRW expansion. Hence, when we carry out the 
flts, we can approximately marginalize over the spectral 
index (since each constant spectral index corresponds to 
a different inflationary model in which the spectral index 
is approximately constant). 

The coefficient A in this approximation of negligible 
back reaction can be related directly to the coupling con- 
stant A by noting that 



This give us 



A = 



dkpki |/3fe| 
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(12) 



2^/b^TV^ VMMU^ 



(13) 
(14) 



where we have used the usual approximation for the pri- 
mordial slow roll inflationary spectrum [5]. This means 
that regardless of the exact nature of the inflationary 
scenario, for any flxed inflationary spectrum (5//(fc)|A=o 



without the back reaction, we have the particle produc- 
tion giving us a bump of the form Eq. (9) with the pa- 
rameter A expressed in terms of the coupUng constant 
through Eq. (14). Given that the CMB normalization 
requires SH{k)\\=o ^ 10~^, we have 



A- 1.3iVA^/^ 



(15) 



Hence, for A ^ 0(0. 1), both A < 1 and AA^ < 1 are 
possible, satisfying perturbativity. As will be seen in the 
sections below, our best fits will indeed give A ^ 0(0. 1). 
When TV and A are sufficiently large that the back re- 
action becomes important (when the naively computed 
A ~ 0(3)),^ the back reaction reduces [22] the actual 
amplitude of the peak relative to the perturbative value 
given by Eq. (14). 

Before we conclude this section, we would like to ex- 
plicitly state further caveats to using Eq. (2) in comput- 
ing the primordial perturbation spectrum: 

1 . Instead of the usual long wavelength gauge invari- 
ant curvature perturbation variable implicit in us- 
ing Eq. (2), the gauge invariant gravitational poten- 
tial must be recomputed incorporating the coupling 
of Eq. (3). 

2. Even with the usual gauge invariant curvature per- 
turbation formalism, Eq. (2) neglects the change 
in pressure that occurs due to the particle produc- 
tion. Even considering the adiabatic perturbation 
component, this would change the denominator of 
Eq. (2). 

3. The numerator of Eq. (2) also obtains a contribu- 
tion from the particle production which we are ne- 
glecting. 

4. As pointed out by [23] , perhaps the most significant 
effect that has been neglected is the nonadiabatic 
pressure change. 

A more detailed investigation of these neglected effects 
will be deferred to a future publication. 



III. MATTER POWER SPECTRUM 

It is straight forward to determine the matter power 
spectrum to compare with that deduced from large-scale 
structure surveys [17,18] and the Lyman-a forest [19]. 
To convert the amplitude of the perturbation as wave 
number k enters the horizon, (5/f (/c), to the present-day 
power spectrum, P{k), which describes the amplitude of 
the fluctuation at a fixed time, one must make use of 



a transfer function, T{k) [29] which is easily computed 
using the code CMBFAST [28] for various sets of cos- 
mological parameters (e.g. Vt, Hq, A, ils)- An adequate 
approximate expression for the structure power spectrum 
is 



7,3 



aHn 



T\k)5l{k) . 



(16) 



^The back reaction becomes very strong as yl ^ 3e ~ 8. 



This expression is only valid in the linear regime, which 
in comoving wave number is up to approximately A:^0.2 h 
Mpc~^ and therefore adequate for our purposes. How- 
ever, we also correct for the nonlinear evolution of the 
power spectrum [30]. 



IV. CMB POWER SPECTRUM 

Features in the primordial power spectrum will also ap- 
pear in the observed CMB temperature fluctuations. The 
connection between the resonant particle creation and 
CMB temperature fluctuations is straightforward. As 
usual, temperature fluctuations are expanded in spheri- 
cal harmonics, 5T/T = J2i J2m "'imYim{d, <P) {2 < I < oo 
and —l<m<l). The anisotropics are then described 
by the angular power spectrum, Ci — (|a;m|^), as a func- 
tion of multipole number I. One then merely requires the 
conversion from perturbation spectrum P{k) to angular 
power spectrum Ci. This is also easily accomplished us- 
ing the code CMBFAST [28]. As input to CMBFAST we 
adopt the usual power law primordial spectrum plus the 
perturbation due to resonant particle production. For 
speed, CMBFAST does not compute all Ci, but uses a 
spline flt to interpolate. We have checked the stability 
of the CMBFAST results when adding such features to 
the power spectrum by increasing the number of Ci ex- 
plicitly computed. The results were convergent even for 
the default spline fits of CMBFAST. When converting to 
the angular power spectrum, the amplitude of the nar- 
row particle creation feature in SH{k) is spread over many 
values of /. Hence, the particle creation feature looks like 
a broad peak which is easily accommodated even when 
implementing spline fits. 



SZ effect 

It has been proposed [13,14] that the observed spec- 
trum at high I is explained by the SZ effect. However, 
the amplitude of the SZ contribution to the power spec- 
trum is very sensitive to the parameter as describing rms 
mass fluctuation within a flducial 8 h^^ Mpc sphere. The 
amplitude of the expected SZ peak scales as ct|. 

Explaining the excess power in the observed CMB 
spectrum requires that the mass density fluctuation am- 
plitude parameter, as, be slightly above unity [14]. How- 
ever, a variety of independent measures including weak 



lensing [31], galaxy velocity fields [32], galaxy clusters 
at high redshift [33], X-ray emitting clusters [34], and 
the independent value obtained by the WMAP fit [6] to 
the lower multipole data all favor a mean value in the 
range tJs ~ 0.7 — 0.9 (see Ref. [13] for a recent review). 
Note that even a slight reduction of the amplitude pa- 
rameter by 10% is sufficient to reduce the magnitude of 
distortion by more than a factor of two. For the present 
purposes we adopt a conservative value of erg = 0.9 ± 0.1 
as a reasonable prior distribution based upon the various 
independent measures of erg . 

In what follows we fit the amplitude of the SZ contri- 
bution using the SZ power spectrum calculated in [14] 
based upon the analytic halo formalism [35]. This ana- 
lytic form has been shown [13,14] to adequately represent 
the power spectra deduced from numerical simulations. 



V. RESULTS 

We have made a multi-dimensional Markov Chain 
Monte-Carlo analysis [36,37] of the mass power spectrum 
based upon the combined 2dF, SDSS, and Lyman-a 
data in models with and without resonant particle cre- 
ation during inflation to alter the primordial fluctua- 
tion spectrum. We also have independently analyzed the 
CMB using the combined WMAP, CBI, ACBAR, and 
VSA data. In addition, we made a total analysis based 
upon the combined mass and CMB power spectra. 

For simplicity and speed in the present study we 
only marginalized over the five parameters which do 
not alter the matter or CMB transfer functions. 
Hence, the set of free parameters in the analysis is 
(us. As, log (fc*). A, Asz), where Us is the spectral index. 
As is the overall amplitude of the primordial power spec- 
trum, and Asz is the overall amplitude of the SZ effect 
which we relate to an effective cf ■^ parameter for the SZ 
effect. The true ug parameter, however, is determined 
by As once the power spectrum and transfer functions 
are fixed. As usual, both n^ and As are normalized at 
k = 0.05 Mpc^^. As noted above, we adopt a conserva- 
tive prior of erf ^ = 0.9 ± 0.1 as opposed to the best fit 
combined IFM^P analysis value of erg = 0.84 ±0.04. For 
this illustrative study all other parameters for this analy- 
sis were fixed at the optimum W^AMP parameters [6], i.e. 

{h,nbh^,nmh^,nA,T) = (0.71,0.0224,0.135,0.75,0.17). 

One expects that different results would arise by optimiz- 
ing over other parameters, e.g. h, Q,bh? , etc. However, for 
the present exploratory study, we only wish to identify 
angular scales which might be of interest for future study. 
For this purpose, a simple five parameter marginaliza- 
tion is adequate since such parameters do not produce 
the kind of spectrum bump of interest here. 



A. Matter power spectrum fit 

Contours of constant goodness of fit in the A vs. fc, 
plane consistent with the matter power spectrum con- 
straint are shown on Figure 3. These fits are based 
upon the 2dF, SDSS, and Lyman-a data described above. 
Window functions required for fits to these data sets are 
given in approximate analytic form in [18] for the 2dF 
data, and in [17] for the SDSS data. As described in [38], 
a simple x^ is used to fit for the Lyman-a data, because 
the full correlation matrix is not available. 

The most recent SDSS matter power spectrum extends 
into the nonlinear regime, but is probably not reliable 
[3] for the last band widths above k = 0.2 h Mpc^^. 
Hence, we have omitted the last points of the SDSS for 
k>0.2 h Mpc-i and for k > 0.15 h Mpc^^ for the 2dF 
power spectrum as recommended [3] . For the SDSS and 
2dF power spectra the unknown bias factors (multiplier 
to get from the matter power spectrum to the galaxy 
power spectrum) were analytically marginalized accord- 
ing to the method described in the appendix of Ref. [37]. 

The contours on Figure 3 identify two possible regions 
which could be consistent with resonant particle produc- 
tion. The feature with the largest amplitude occurs at 
k^ = 0.41±0.05/iMpc~i and with yl = 0.23±0.08. These 
parameters correspond to tti w 2.2 M^; and A w 0.6 for a 
single fermion species. 

The weaker feature has a minimum at k^ — 0.08 ± 
0.01 h Mpc-i and an amplitude A = 0.12 ± 0.04. This 
feature is more or less consistent with a similar feature 
seen in the CMB, but at a lower amplitude and slightly 
larger angular scales. This feature would correspond to 
m « 1.8 Mpi and A « 0.5. However, it occurs at the 
interface between the SDSS and Lyman-a power spec- 
tra, and hence, may be an artifact of the matching of 
the data sets rather than a real bump in the spectrum. 
Obviously, this is another feature which warrants further 
careful scrutiny. 
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FIG. 3. Constrains on parameters A and k, from the fit to 
the matter power spectrum alone. Contours show I and 2 a 
confidence limits. The horizontal axis indicates log(fe,) where 
fc, is in units of {h Mpc"^). 

Figure 1 shows the effect of the stronger feature on the 
matter power spectrum compared with data from 2dF, 
SDSS, and Lyman-a. Clearly the feature in the mat- 
ter power spectrum beginning near k « 0.6 h Mpc "^ is 
important in this regard and warrants careful scrutiny. 
(Note that the peak of the primordial spectrum bump 
generically occurs at fc w cxp(l/3)fc, with our parame- 
terization.) 



B. CMB Fit 

Cosmological parameters with and without resonant 
particle creation were obtained from the combined 
WMAP, CBI, ACBAR, and VSA data using a Markov 
chain Monte-Carlo analysis [36] as described above. Fig- 
ure 4 shows the 1 and 2ct contours in the fc, vs. A plane 
from this analysis. The CMB power spectrum is best fit 
for A = 0.7± 0.2 and fc, = 0.18 ± 0.02 h Mpc^^. 

Figure 2 shows our best fit to the CMB data for models 
with and without resonant particle creation. Also shown 
for comparison in the lower expanded figure is the SZ 
contribution for erg = 0.9. We find that an excellent fit to 
the ACBAR and CBI observed CMB power spectra can 
be achieved in this way. Indeed, those CMB observations 
favor this interpretation over the conventional SZ plus 
kCDM models at the level of ha. However, the BIMA 
data set at I k, 5000—7000 is a crucial test of this possible 
interpretation. If the power spectrum is as high at these 1- 
values as the BIMA data imply, then an SZ interpretation 
with a large cg is favored. 
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FIC. 4. Constrains on parameters A and fc, from the CMB 
power spectrum. Contours show 1 and 2(7 limits. The hori- 
zontal axis indicates log(fc,) where fc* is in units of [h Mpc^^). 



C. Combined Analysis 

Contours of A and fc.^ consistent with the combined 
CMB and matter power-spectrum constraints are shown 
on Figure 5. Here it is apparent that a single fea- 
ture which begins in the matter power spectrum at 
A:* = 0.17 ± 0.04 /iMpc^^ and has an amplitude of 
A K, 0.35 ± 0.10 is most prominent. Figure 6 illustrates 
the associated optimum fits to the matter power spec- 
trum (upper curve) and the CMB (lower curve). 

As remarked above, even though a feature with fc, ~ 
0.2 appears in both the matter and CMB power spectra, 
it is not unambiguously attributable to resonant parti- 
cle production. In the matter power spectrum it could 
be an artifact of the uncertainty in joining of the galaxy 
and Lyman-a data sets near k « 0.1, while in the CMB 
its significance could be confused by the SZ effect. Never- 
theless, the feature with fc* = 0.4 ft,Mpc^^ and A k, 0.2 
which is clearly seen in the matter power spectrum of 
Figure 1 is still apparent in the combined analysis. The 
diminished significance in the combined analysis is to be 
expected since this bump is entirely due to the Lyman-a 
data and is not detectable in the CMB. 

Regarding other parameters in the analysis, it is not 
surprising that we deduce values very near to the WMAP 
[6] optimum results, e.g. Ug — 0.96 ± 0.01, a^ — 0.88 ± 
0.08. 
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FIG. 5. Constrains on parameters A and k^ from a fit to 
the combined CMB and matter power spectra. Contours show 
1 and 2a hmits. The horizontal axis indicates log(/i;*) where 
k^ is in units of {h Mpc~^). 
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FIG. 6. Optimum fit to the combined CMB -I- matter 
power spectra. Upper figure sfiows the observed galaxy clus- 
ter function from SDSS [17], 2dF [18], and Lyman-a [19] data. 
The solid line shows the spectrum implied from the optimum 
model with resonant particle creation during inflation. The 
lower figure shows an expanded view of the ACBAR, CBI, 
VSA, and BIMA data in the range of I = 1000 - 7000. The 
dashed line shows the standard WMAP result without the SZ 
effect included. The solid line shows the optimum fit with 
resonant particle production and the SZ contribution. The 
dOt-dot-dashed line shows the resonant particle creation com- 
ponent without the SZ effect. 



VI. CONCLUSION 

We have analyzed the mass and CMB power spectra 
in the context of a model for particle creation during 
inflation. We find marginal evidence for excess power 
in both the CMB and mass fluctuation spectra consis- 
tent with this hypothesis. The combined CMB -I- mat- 
ter power spectrum would imply an optimum feature at 
K = 0.17 ± 0.04 /iMpc-i and A « 0.35 ± 0.10. How- 
ever, given that the CMB results may be attributed to 
the SZ effect and that the matter power spectrum is very 
uncertain at this scale, we propose that the most likely 
possibility could be the feature in the matter power spec- 
trum deduced from the Lyman-a absorption structures 
at a scale of /c ~ 0.4 h Mpc~^ and A « 0.2. Either of 
these features would correspond to the resonant creation 
of a particle with m w 1 — 2 M^; and a Yukawa coupling 
constant between the fermion species and the inflaton of 
A ~ 0.5 — 1.0 for a single fermion species. 

Obviously there is a need for more precise determina- 
tions of the primordial matter power spectrum on the 
scale of 0.1 to 10 Mpc. Indeed, the Lyman-a forest con- 
straints and method are currently under revision and 
new results are expected shortly which may alter the 
conclusions presented here. There is also a need for 
more precise determinations of the SZ contribution to 
the CMB on the scale of I « 2000 - 10000. We further 
note that the present analysis has neglected the possibil- 
ity of nonadiabatic effects. Based upon previous stud- 
ies [39,40] of nonadiabatic isocurvature fluctuations on 
the matter power spectrum and CMB, we would expect 
that the introduction of non-gaussian isocurvature fluc- 
tuations could add excess power on small angular scales 
and hence could decrease the significance of the features 
identified here. They would not, however, naturally pro- 
duce the kind of bump of interest here, and hence would 
not easily explain such features away. 

In spite of these caveats, we conclude that if the present 
analysis is correct, this may be one of the first hints at ob- 
servational evidence of new particle physics at the Planck 
scale. Indeed, one expects a plethora of particles at the 
Planck scale, particularly in the context of string theory, 
the leading candidate for a consistent theory of quantum 



gravity. Perhaps, the presently observed power spectra 
contain the first suggestion that such particles may have 
not only existed in the early universe, but coupled to the 
inflaton field, and thereby, left a relic signature of their 
existence in the primordial power spectrum. 
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